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r^ • Abstract 

Fidelity is known to increase through any Kraus map: the fidehty between two 

G I density matrices is less than the fidelity between their images via a Kraus map. 

^ ■ We prove here that, in average, fidelity is also increasing for discrete-time quantum 

O^i filters attached to an arbitrary Kraus map: fidelity between the density matrix of 

the underlying Markov chain and the density matrix of the associated quantum filter 

cn ■ is a sub-martingale. This result is not restricted to pure states. It also holds true for 

j^ ■ mixed states. 

^ '. 1 Introduction 
in 

Q ■ Quantum filtering was developed, in its modern theoretical form, by Belavkin during the 

1980's ['.] (see also [2] for a recent tutorial introduction). Quantum filtering proposes a 
theory for statistical inference strongly inspired from quantum optical systems. These sys- 

k> ! tems are described by continuous-time quantum stochastic differential equations (stochastic 

^ [ master equations). The stability of the obtained quantum filters, i.e, the fact that quan- 

tum filters are asymptotically independent of their initial states, has been investigated in 
several papers. In [9] sufficient convergence conditions are established: they are related 
to observability issues.^ In [i] convergence of the estimate to the physical state is ad- 
dressed for a generic Hamiltonian and measurement operator: the analysis relies on the 
fact that, for pure states, the fidelity between the real state and its estimate is proved to 
be a sub-martingale. 

Quantum filtering is not restricted to continuous time. It can also be used for discrete- 
time systems as in [5] where a quantum filter is used inside a state-feedback scheme. This 
note proves that fidelity between the estimate and real state is always a sub-martingale 
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^See [1(J] for general connections between the stability of nonlinear filters and mathematical systems 
theory concepts such as stability, observability and detectability. 



for such discrete-time quantum filters. This does not mean that quantum filters are al- 
ways convergent. Nevertheless, this sub- martingale combined with stochastic invariance 
arguments could be useful to prove convergence in specific situations. 

Fidelity F{a, p) between two quantum states described by density matrices p and a 
coincides with their Frobenius inner product Tr {per) when at least one of two density 
matrices p or cr is a projector of rank one, i.e. a pure state. ^ The fact that fidelity 
is a sub martingale for pure states is a known fact. For continuous-time systems, the 
convergence investigations proposed in [1] relies on this fact. For discrete-time systems, 
theorem 3 in [7], based only on Cauchy-Schwartz inequalities, proves that the Frobenius 
inner product between the estimate and real state is always a sub-martingale whatever the 
purity of the estimate and real state are. Consequently, if we assume that the real state 
remains pure, then, fidelity coincides with Frobenius inner product and thus is known to 
be a sub- martingale. When the real state and its estimate are mixed, fidelity does not 
coincide with Frobenius inner product and the proof that it is a sub- martingale does not 
result directly from published papers, as far as we know. Fidelity between the arbitrary 
states p and a is given by Tr^ ( \f'^fpo\fp) , an expression much more difficult to manipulate 
than Tr (per). The contribution of this note (theorems 1 and 2 relying on inequalities (4) 
and (5), respectively) consists in proving that, in the discrete-time case and for arbitrary 
purity, the fidelity between the estimate and real state is sub-martingale. 

In section 2, we recall the Kraus map associated to any quantum channel, the associated 
Markov chain (quantum Monte-Carlo trajectories) and highlight the basic inequalities (4) 
and (5) underlying theorems 1 and 2 proved in the remaining sections. In section 3, we 
consider the quantum filter attached to a quantum trajectory and prove theorem 1: fidelity 
is a sub-martingale. In section 4 we extend this result to a family of Markov chains attached 
to the same Kraus map. In section 5, we conclude by noting that, contrarily to fidelity, 
trace distance and relative entropy are not always super-martingales. 

The author thanks Mazyar Mirrahimi and Ramon Van Handel for fruitful comments 
and interesting references. 

2 Kraus maps and quantum Markov chains 

Take the Hilbert space S = C" of dimension n > and consider a quantum channel 
described by the Kraus map (see [8], chapter 8 or [6], chapter 4) 

m 

/C(p) = J]M,pMt (1) 

where 



^There are two closely related definitions of the fidelity. The "mathematical" definition used, e.g., in [8] 
corresponds to the square root of the "physical" definition we have adopted in this paper for simplicity 
sakes. 



• p is the density matrix describing the input quantum state, /C(p) being then the 
output quantum state; p G C"^" is a density matrix, i.e., an Hermitian matrix semi- 
positive definite and of trace one; 

• for each /i G {1, . . . , m}, M^ G C"^"/{0}, and ^^ MtM^ = /. 

To this quantum channel is associated the following discrete-time Markov chain (quantum 
Monte-Carlo trajectories, see, e.g., [6]): 

Pk+i = M^,{pk) (2) 

where 

• pk is the quantum state at sampling time t^ and k the sampling index (t^ < t^+i). 

• /ifc G {1, . . . , m} is a random variable; pk = fi' with probability p^{pk) = Tr (M^pfc-^/1) ■ 

• MJp) = —r^ — TxM.pM}, = -^M.pMl . 

Quantum probabilities derive from density matrices which are dual (via the Frobenius 
Hermitian product Tr (^AB'^), A, B, n x n matrices with complex entries) to the set of 
Hermitian n x n matrices (space of observables). The Kraus map (1) leads to a gener- 
alized observable (or effect valued measure) V H- J2u£V ^l^t^ where V is any subset of 
{1, . . . , m}. The map (V, p) H- XlueP ^tJ-P-^l ^^ known as Davies instrument (see chapter 
4 of [3]). It formalizes a measurement notion relating the state (the density matrix p), mea- 
sure properties (here V living in the a-algebra generated by the singletons of the discrete 
set {1, . . . , m}) and an output state conditional on the observed values (here J^a^r ^fj^P^l 
up to a normalization to ensure a unit trace). The non-linear operators Ai^ just define the 
one-step transition mechanism for the above quantum Markov chain. Such setting can be 
seen as an appropriate analogue of what one encounters in classical discrete-time filtering 
problems. 

Fidelity is one of several ways to measure difference between density matrices and is 
specific to quantum theory. The Kraus map tends to increase fidelity F (see [8], theorem 
9.6, page 414): for all density matrices p and a, one has 

Tr^ (VT^^^i^ifW^) = F{lC{a),]C{p)) > F{a,p) = Tr^ (jV^pV^) (3) 

where, for any Hermitian semi-positive matrix A = UAW, U unitary matrix and A = 
diag{A/}z6{i,...,n}, VA = UVAU^ with a/A = diag{^/X'l}l(z{l^,„^n}■ 

The conditional expectation of pk+i knowing pk is given by the Kraus map: 

E{pk+i/pk) = IC{pk). 



This result from the trivial identity ^"^^i Tr i^M^pM'l^ tTajHu^' ~ ^{p)- ^^ section 3, we 
show during the proof of theorem (1) the following inequality 

y Tr (M.pMl) F ( -JkzMl^^ f^f^l ] > pia, p) (4) 

for any density matrices p and a. The left-hand side is related to a conditional expectation. 
Inequahty (4), attached to the probabilistic mapping (2), can be seen as the stochastic 
counter-part of inequality (3) attached to the deterministic mapping (1). When for some 
/i, Tr (M^crMt) = with Tr [M^pMl) > 0, the /i-term in the sum (4) is not defined. This 

is not problematic, since in this case, if we replace — . '"'^ ^t\ by — r-^^ — ^ where f is any 

density matrix such that Tr {^M^^M^\ > 0, this term is then well defined (in a multi- valued 
way) and inequality (4) remains satisfied for any such ^. 

During the proof of theorem (8), we extend this inequality to any partition of {1, ... , m} 
into p > 1 sub-sets Vy. 

E T. ( E M,,Mi\ F ( ,|;-X^., . ^isf^) > n. .) (5) 

3 The standard case. 

Take a realization of the Markov chain associated to the Kraus map /C. Assume that 
we detect, for each fc, the jump /i^ but that we do not know the initial state po- The 
objective is to propose at sampling /c, an estimation p^ of p^ based on the past detections 
Po, • • • ,Pfc-i- The simplest method consists in starting from an initial estimation po and 
at each sampling step to jump according to the detection. This leads to the following 
estimation scheme known as a quantum filter (see, e.g., [11]): 

Pk+i = M^,,{pk) (6) 

with Pfj,{pk) = Tr (M^pkM^) as probability of pk = /^- Notice that when Tr {M^^p^M^^) = 
0, J^^^{pk) is not defined and should be replaced by M.fi^.{C) where ^ is any density matrix 
such that Tr {M^^^M^A > (take, e.g., ^ = ^Id)- The theorem here below could be 
useful to investigate the convergence of pk towards pk as k increases. 

Theorem 1. Consider the Markov chain of state {pk,Pk) satisfying (2) and (6). Then 
F{pk,pk) is a sub-martmgale: E {F{pk+i, pk+i) / {pk, Pk)) > F{pk,pk)- 

The proof proposed here below deals with the general case when both pk and pk can be 
mixed states. It relies on arguments similar to those used for the proof of theorem 9.6 in [8]. 
They are based on a kind of lifting process, usual in differential geometry. They consist in 
introducing additional variables and embedding the Hilbert space 5* into the tensor product 



with a copy Q oi S and with the environment-model space E. In the "augmented space" 
S ^ Q ^ E, the quantum dynamics (2) read as an unitary transformation followed by an 
orthogonal projection. At sampling step k, the estimated and real states are lifted to the 
"augmented space" via a standard procedure known as purification. These purifications are 
then subject to an unitary transformation followed by the same projection. Partial traces 
over Q ® E provide then the estimated and real states at sampling time k + 1. At each 
sampling step A;, the purifications of the estimated and real states are chosen according 
to Ulhmann's theorem: the square of modulus of their Hermitian product is maximum 
and coincides thus with the fidelity. The remaining part of the proof is based on simple 
manipulations of Hermitian products and Cauchy-Schwartz inequalities. 

Proof. The density matrices p and p are operators from 5 = C" to S. Take a copy Q = C" 

2 

of S and consider the composite system living on S* ® Q = C" . Then p and p correspond 
to partial traces versus Q of projectors \ip){ip\ and |^)(^| associated to pure states \%Ij) and 
\^) eS®Q: 

p = TrQ(|z^)(^|), p = TrQ(|^)(^|) 

\ip) and \il)) are called purifications of p and p. They are not unique but one can always 
choose them such that F{p,p) = |(?/;|?/^)p (Uhlmann's theorem, see [8], theorem 9.4 page 
410). Denote by \ipk) and \ipk) such purifications of p^ and pk satisfying 

Fipk,pk) = \{i'km\'. 

We have 

m 

E{F{pk+i,Pk+i)/{pk,Pk)) = ^P;.(pfc)F(A^^(pfc),A^^(pfc)). 

The matrices A4fj,{pk) and J^^{pk) are also density matrices. 

We will consider now the space S ^ Q ^ E where E = C" is the Hilbert space of the 
environment appearing in the system-environment model of the Kraus map (1) (see [6], 
chapter 4 entitled "The environment is watching"). Thus exist |eo) a unitary element of 
E, an ortho-normal basis (|p))^e{i,...,m} of E and a unitary transformation U (not unique) 
on S ^ E such that, for all |0) G S", 

m 

For each p, denote by P^ the orthogonal projector onto the subspace S (8) (C|p)). Then 
P^U (10) (g> |eo)) = (Af^|0)) ® |p) and XIu-^m ~ ^- -^^^ ^^y density matrix p on S, 

P,U (p ® |eo)(eo|) U^P, = M^pMl ® |p)(p| 

and thus 

Tr^ {P^U (p ® |eo)(eo|) U^P,,) = M^pMl 



The unitary transformation U oi S ^ E can be extended to S^Q®E = S^E®Qhj 
setting V = U ^ I {I is identity on Q). Since \ipk) ^ S ^ Q is a. purification of pk, we have 

M^PkMl = Tte^q {P,V{\^k){i^k\ ® \eo){eo\)V^P,) . 

Set |0fc) = lipk) ® |eo) e S ^Q E and \xk) = V'l^fc). Using P^ = P^, we have 



p,{pk) = Tr (M^PfcMt) = {<Pk\V^P,V\<Pk) = \\P,\Xk)\\ 
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Thus, for each u, the state \xka) = / ^ PulXk) is a purification of MJpk)- 



M^,{pk) = Trq^E ilXkf,) {Xk^,\) ■ 

eo) and \xk) = V\(t>k)- F( 
also a purification of Ai^{pk)- By Uhlmann's theorem 



Similarly set |0fc) = \^k) ® |eo) and \xk) = V\(t>k)- For each p, \xkfi) = -^=^Pf,\Xk) is 
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Thus we have 

m 
E{F{pk+l,pk+l)/{pk,Pk)) > ^P,^{pk) liXkf^lXkf,] 

Since V is unitary, 

\{xk\xk)\' = \{k\<Pk)\' = mmi' = F{p,,p,). 

Let us show that Yl'i^^iPpiPk) l(Xfc/.|XfcM)P ^ \{Xk\Xk)\^- We have 

Pfiipk) \{Xkt,\Xkt,)\'^ = \{Xkt,\P^Xk)\'^ = \{Xk^\Xk)\'^, 

thus it is enough to prove that J2^=i \{Xkfi\Xk)\'^ > |(Xfc|Xfc)P- Denote hjRcS^Q^E 
the vector space spanned by the ortho-normal basis (|XfcM))ue|i ml ^^^ ^^ ^^^ projector 
on R. Since 



\xk) = Xl^^l^fc) = XI V^^'^^'^)!^'^^^ 

IXfc) belongs to R and thus |(Xfc|XA;)P = \{Xk\P\Xk))\^- We conclude by Cauchy-Schwartz 
inequality 



\{xk\xk)\' = \{xk\P\xk))\' < \\xkf\\P\xk)f = \\P\xk)\\' = 5Z Kxfc^lx;^-) 

p=l 



|2 



n 



4 The coarse-grained case. 

Let us consider another Markov chain associated to a partition of {1, . . . ,m} into p > 1 
sub-sets Vu (aggregation of several quantum jumps via "partial Kraus maps"): 



Pk+l 



Tr(E,e.,, M.P.M., 



TY E ^'^P'^K (7) 



where Vk = i^ with probability Tr ( ^^g-p^ Mfj_pkMj^ j . It is clear that E (pfc+i/pfc) = ^{pk)'- 

the Markov chains (2) and (7) are attached to the same Kraus map (1). Consider the 
associated quantum filter 

P^+' = Tf¥ 'm ' mA I E M.PkMl (8) 

where the jump index Uk coincides with the jump index z/^ in (7). Then we have the 
following theorem. 

Theorem 2. Consider the Markov chain of state {pk,Pk) satisfying (7) and (8). Then 
F{pk,pk) is a sub-martingale: E {F{pk+i, pk+i) / {pk, Pk)) > F{pk,pk)- 

Proof. It is similar to the proof of theorem 1. We will just point out here the main changes 
using the same notations. We start from 

p 
E{F{pk+i,pk+i)/{pk,Pk)) = '^Pu{pk)F{M^{pk),Mu{pk))- 

v=l 

where we have set 

P.(p) = Tr ( 5^ M.pMtj , ^.(p) = ^($^M,pMt). 

With Pi, the orthogonal projector on S^Q®span{\p), p G V^} and My{p) = J2uev^ ^nP^l^ 
we have 

J2 M^pkMl = TiQ^E (PuV {\ilJk){M ® |eo)(eo|) V^P,'^ 

and 

P,{pk) = Tr ( J] M^pkMU = {<Pk\V^P.V\<Pk) = \\Pu\Xk) 

For each i/, the state \xkv) = , ^ PulXk) is a purification of Aiu{pk)'- 

Muipk) = Trgcgij i\Xku){Xk,y\) ■ 



Similarly \xki) = , "^ . PvYkk) is also a purification oi M.y{pk). By Uhlmann's theorem, 

Thus we have 

p 

E(F(pfc+i,Pfc+i)/(pfc,Pfc)) >'^Pu{pk) \{Xk^,\Xk^J,)\'^■ 

Let us show that YZ=iPy{pk) \{Xku\Xku)? > |(Xfc|Xfe)P = F{pk,pk)- We have 
Puipk) \{Xku\Xku)\'^ = \{Xku\PuXk)\'^ = \{Xku\Xk)\'^, 

thus it is enough to prove that Yl^=i \{Xku\Xk)\'^ > |(Xfe|Xfc)P- Denote hjRcS<^Q®E 

the vector space spanned by the ortho- normal basis ( \xku) ) and by P the projector 

V / i/e{i,...,p} 

on R. Since 

\Xk) = X^^^IXfc) = X^ VpjJkilXki') 

u=l v=l 

\Xk) belongs to R and thus |(Xfc|XA:)P = \{Xk\P\Xk))\'^ ■ We conclude by Cauchy-Schwartz 
inequality 

p 

\{Xk\Xk)? = \{Xk\P\Xk))\' < \\Xkf\\P\Xk)f = \\P\Xk)\\' = J2\{Xku\Xk)\'- 

n 



5 Concluding remarks 

Theorems 1 and 2 are still valid if the Kraus operators M^ depend on k. In particular, 
F{pk,Pk) remains a sub-martingale even if the Kraus operators depend on p^, i.e., in case of 
feedback. A natural extension will be to prove that fidelity yields also to a sub-martingale in 
the continuous-time case. The Markov chain (2) is then replaced by a stochastic differential 
equation for p driven by a Poisson or a Wiener process. 

Kraus maps are contractions for the trace distance (see [8], theorem 9.2, page 406): for 
all density matrices a, p, one has Tr(|/C((T) — /C(p)|) < Tr(|(T — p|). When a and p are 



pure states (projectors of rank one), D{a, p) = a/1 — F{a,p). Consequently inequality (4) 
yields to 

y IV (M,pMl) D I ^lii^, _Jft£Ml ] < D{a, p) 



^l=l 



for any pure states a and p (use the fact that [0, 1] 3 x H- ^Jl — x is decreasing and 
concave). It is then tempting to conjecture that the above inequality holds true for any 
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mixed states p and a, i.e., that D{pk,pk) = Tr(|pfc — pfc|) is a super-martingale of the 
Markov chain defined by (2) and (6). Unfortunately, this is not true in general as shown 
by the following counter-example. Take n = 3, m = 2 and 




a 




0\ /O 

M, = ( ^ , M,= [0 ^ 
0/ \0 1, 



Then E,t,Tr(M,.pAf.)fl(,j^.5Mg|^) = f > 1 = D(., p).' The same 

counter example shows that the relative entropy Tr(pyfclogpfc) — Tr (p^,. logp^) is not, in 
general, a super-martingale. 
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